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Abstract 

The  proton  exchange  membrane  (PEM)  fuel  cell  is  being  investigated  as  an  alternate  power  source  for  various  applications  like  transportation 
and  emergency  power  supplies.  The  paper  presents  a  novel  circuit  model  for  a  PEM  fuel  cell  that  can  be  used  to  design  and  analyze  fuel 
cell  power  systems.  The  PSPICE-based  model  uses  bipolar  junction  transistors  (BJTs)  and  LC  elements  available  in  the  PSPICE  library  with 
some  modification.  The  model  includes  the  phenomena  like  activation  polarization,  ohmic  polarization,  and  mass  transport  effect  present  in 
a  PEM  fuel  cell.  The  static  and  dynamic  characteristics  obtained  through  simulation  are  compared  with  experimental  results  obtained  on  a 
commercial  fuel  cell  module. 
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1.  Introduction 

Fuel  cells  of  various  types  are  considered  as  alternatives  to 
fossil  energy  mostly  for  reasons  of  pollution  and  efficiency. 
The  proton  exchange  membrane  fuel  cell  (PEMFC)  has  been 
considered  as  a  promising  kind  of  fuel  cell  during  the  last  10 
years  because  of  its  low  working  temperature,  compactness, 
and  easy  and  safe  operational  modes.  The  proton  exchange 
membrane  (PEM)  fuel  cell  is  very  simple  and  uses  a  polymer 
(membrane)  as  the  solid  electrolyte  and  a  platinum  catalyst. 
The  hydrogen  from  a  pressurized  cylinder  enters  the  anode 
of  the  fuel  cell  and  the  oxygen  (from  air)  enters  the  cathode. 
Protons  and  electrons  are  separated  from  hydrogen  on  the 
anode  side.  In  a  basic  PEM  cell,  the  protons  are  transported 
to  the  cathode  side  through  the  polymer  and  the  electrons 
are  conducted  through  the  load  outside  the  electrode  [1].  A 
fuel  cell  stack  is  composed  of  several  fuel  cells  connected  in 
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series  separated  by  bipolar  plates  [2]  and  provides  fairly  large 
power  at  higher  voltage  and  current  levels. 

There  is  a  need  to  model  the  PEMFC  for  optimizing  its 
performance  and  also  for  developing  fuel  cell  power  convert¬ 
ers  for  various  applications.  Almost  all  the  models  proposed 
for  the  PEMFC  consist  of  mathematical  equations  and  are  not 
of  much  use  in  power  converter/system  simulation  and  anal¬ 
ysis  [1-4].  Other  models  of  PEMFC  use  Matlab-Simulink 
[5]  and  PSPICE  [6],  but  they  are  still  mathematical  in  nature. 
The  models  include  several  chemical  phenomena  present  in 
the  fuel  cell  and  hence  are  complex.  Some  of  the  physical 
variables  like  pressure  and  hydrogen  input  are  constrained  in 
a  commercial  fuel  cell  module  and  this  makes  the  fuel  cell  op¬ 
eration  simpler.  This  also  allows  the  use  of  a  simpler  electric 
circuit  model  useful  to  a  power  electronics  designer. 

This  paper  presents  a  novel  circuit  model  for  a  PEMFC, 
which  is  simple  and  at  the  same  time  includes  all  the  important 
characteristics  of  a  fuel  cell  stack.  The  model  uses  the  non¬ 
linearity  of  a  junction  diode  and  the  current  control  feature 
of  bipolar  junction  transistors  (BJTs).  In  the  proposed  model, 
a  diode  is  used  to  model  both  the  activation  losses  and  the 
ohmic  losses  in  a  PEMFC,  while  two  BJTs  are  used  to  model 
the  mass  transport  losses  [7,8].  The  equations  governing  the 
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different  polarization  effects  in  a  fuel  cell  are  related  to  the 
equations  of  the  electrical  circuit  elements  making  up  the 
model.  The  model  is  justified  through  experiments  performed 
on  a  commercial  PEMFC  stack  for  both  static  and  dynamic 
responses. 


2.  Polarization  characteristics  of  a  PEM  fuel  cell 

Proton  exchange  membrane  fuel  cells  combine  hydrogen 
and  oxygen  over  a  platinum  catalyst  to  produce  electrochem¬ 
ical  energy  with  water  as  the  byproduct.  Fig.  1  shows  the 
V-I  characteristic  of  a  typical  single  cell  operating  at  room 
temperature  and  normal  air  pressure  [9] .  The  variation  of  the 
individual  cell  voltage  is  found  from  the  maximum  cell  volt¬ 
age  (or  EMF)  and  the  various  voltage  drops  (losses).  Multiple 
factors  contribute  to  the  irreversible  losses  (voltage  drop)  in 
an  actual  fuel  cell  that  cause  the  cell  voltage  to  be  less  than 
its  ideal  potential  [9].  The  losses,  which  are  also  called  polar¬ 
ization,  originate  primarily  from  three  sources:  (a)  activation 
polarization,  (b)  ohmic  polarization,  and  (c)  concentration 
(mass  transport)  polarization.  Each  of  these  is  associated  with 
a  voltage  drop  and  is  dominant  in  a  particular  region  of  cur¬ 
rent  density  (low,  medium,  or  high).  Fig.  1  shows  the  different 
regions  and  the  corresponding  polarization  effects. 

The  ideal  voltage  is  the  maximum  voltage  that  each  cell 
in  the  stack  can  produce  at  a  given  temperature  with  the  par¬ 
tial  pressure  of  the  reactants  and  products  known.  For  the 
PEMFC,  where  pure  hydrogen  and  air  are  used,  the  ideal 
voltage  can  be  calculated  based  on  Gibbs  free  energy  and 
it  is  equal  to  1.2  V  at  25  °C  and  atmospheric  pressure  for  a 
single  fuel  cell  [9].  A  higher  output  voltage  is  obtained  by 
connecting  several  cells  in  series.  The  area  of  the  cell  decides 
the  output  current. 

2.7.  Activation  polarization 

The  activation  polarization  loss  (dominant  at  low  current 
density)  is  present  when  the  rate  of  the  electrochemical 
reaction  at  an  electrode  surface  is  controlled  by  sluggish 
electrode  kinetics  [9].  Activation  losses  increase  as  the 
current  increases.  The  activation  losses  can  be  obtained  by 
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Tafel  equation  [8]: 

V  =  A\n(j2j  (1) 

where  A  is  the  constant,  V  the  over-voltage,  i  the  current 
density,  and  i0  is  the  current  density  at  which  the  voltage 
begins  to  drop. 

2.2.  Ohmic  polarization  (loss) 

The  ohmic  loss  is  due  to  the  resistance  of  the  polymer 
electrolyte  membrane  to  the  ions  and  the  resistance  of  im¬ 
perfect  electrodes.  The  loss  (voltage  drop)  in  the  fuel  cell  is 
approximately  linear  in  this  region. 

23.  Concentration  polarization  ( mass  transportation 
losses ) 

The  concentration  polarization  relates  to  the  change  in  the 
concentration  of  the  reactants  at  the  surface  of  the  electrodes 
as  the  fuel  (hydrogen)  is  used.  The  concentrations  of  the  fuel 
and  oxidant  are  reduced  at  the  various  points  in  the  fuel  cell 
gas  channels  and  are  less  than  the  concentrations  at  the  inlet 
portion  of  the  stack.  This  loss  becomes  significant  at  higher 
currents  when  the  fuel  and  oxidant  are  used  at  higher  rates 
and  the  concentration  in  the  gas  channel  is  at  a  minimum. 

In  general,  the  mass  transportation  (transfer)  losses  are 
given  by  Eq.  [7] 

V  =  V'  -  rnlmenlm  (2) 

where  Im  =  I  — 1\ ,  and  V(,  and  I\  are  the  coordinates  of  the 
point  where  the  V-I  characteristic  starts  to  deviate  from  being 
linear  (start  of  mass  transport  action),  and  m  and  n  are  mass 
transfer  parameters.  While  I\  is  the  limiting  current  at  which 
the  fuel  is  used  up  at  a  rate  equal  to  its  maximum  supply 
rate,  other  constants  depend  on  the  fuel  cell  and  its  operating 
condition. 


3.  Circuit  model  of  PEM  fuel  cell 

Fig.  2  shows  the  proposed  circuit  model  of  a  commer¬ 
cial  PEMFC  module.  The  complete  model  is  developed  by 
modeling  the  different  operating  regions  using  elements  from 


Q1  Ls  R2  Ds 


Fig.  1.  V-I  characteristic  of  a  single  PEM  fuel  cell. 


Fig.  2.  Proposed  circuit  model  of  PEM  fuel  cell  module. 
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the  PSPICE  simulation  library.  Both  static  and  dynamic  re¬ 
sponses  are  considered  in  the  model. 

3.1.  Static  model 

3.1.1.  Activation  polarization 

In  the  model  proposed,  a  diode  is  used  to  model  the  acti¬ 
vation  polarization  of  the  fuel  cell.  It  is  well  known  that,  in 
a  semiconductor  diode,  the  depletion  region  provides  a  po¬ 
tential  barrier  to  inhibit  the  migration  of  carriers  across  the 
junction,  which  is  similar  to  the  “sluggish  electrode  kinet¬ 
ics”  in  the  fuel  cell.  Actually,  the  similarity  can  be  seen  by 
comparing  Eq.  (1)  and  the  diode  equation: 

//D\  kT 

VD=rcVrln(y),  VT  =  —  (3) 

where  Vd  is  the  voltage  across  the  diode  and  /d  is  the  cur¬ 
rent  through  the  diode,  n  the  empirical  constant  between  1 
and  2,  /s  the  reverse  saturation  current,  k  the  Boltzmann’s 
constant,  T  the  absolute  temperature,  and  q  is  the  electronic 
charge.  Thus,  it  is  reasonable  to  use  a  diode  to  model  the  ac¬ 
tivation  polarization  in  a  fuel  cell.  In  the  simulation  package 
PSPICE,  the  diode  model  can  be  adjusted  to  match  the  V-I 
characteristics  of  a  fuel  cell  stack  by  choosing  suitable  values 
for  the  following  parameters:  /s  (saturation  current),  Rs  (para¬ 
sitic  resistance),  A  (emission  coefficient),  /kf  (high-injection 
“knee”  current),  /sr  (recombination  current  parameter),  and 
/bv  (reverse  breakdown  “knee”  current)  [11].  Among  the  sev¬ 
eral  parameters,  A  affects  the  shape  of  the  characteristic  most 
and  it  has  to  be  chosen  to  match  the  V-I  characteristics  of  the 
diode  and  the  fuel  cell  stack.  Fig.  3  shows  the  V-I  character¬ 
istics  for  different  values  of  A  along  with  the  experimental 
characteristic.  It  is  seen  that  the  characteristic  with  A=  80 
matches  the  experimental  characteristic  closely.  /s  was  se¬ 
lected  as  0.02  A. 

3.1.2.  Ohmic  polarization 

A  resistor  is  used  to  model  the  ohmic  polarization.  Instead 
of  using  a  separate  resistor,  the  “parasitic  resistance  (. Rs )”  in 
the  diode  is  used  in  the  model. 


3.1.3.  Concentration  polarization 

To  model  the  concentration  or  mass  transport  over-voltage, 
a  “current  limiting  circuit”  is  used.  In  Fig.  2,  the  current  limit¬ 
ing  circuit  is  composed  of  two  B  JTs  Q\ and  02  and  the  current 
sensing  resistor  R2.  When,  the  current  through  R2  exceeds  a 
set  limit,  Q2  starts  conducting  and  reduces  the  base  voltage 
of  Q\ .  As  a  result,  the  emitter  voltage  of  Q\  (output  voltage) 
will  decrease  at  an  exponential  rate  similar  to  Eq.  (2).  The 
variation  of  output  voltage  as  a  function  of  load  current  is 
obtained  using  the  steady  equivalent  circuit  shown  in  Fig.  4. 
The  transistors  Q\  and  Q2  are  assumed  to  be  identical  with 
current  gain  /3  and  base-emitter  voltage  Vbe-  The  following 
equations  can  be  written: 


Ib\ 


Iq  ~  lei 

1+0 


IC2  ~  IcS^IqR2^Vt 


(4) 

(5) 


Vo  =  Vs  -  Ri  (/bi  +  I  Cl)  -  VbeI  -  (Io  -  I  Cl  -  Ibi)Ri 


(6) 


Substituting  (4)  and  (5)  in  (6)  and  simplifying,  we  get 


V0  =  Vs  -  I, 


(*2+I 


R\ 


+  P 


—  Vbe  — 


Ri  1- 


1 


+R2 


±1 

P 


fee 


I0R2/Vt 


Assuming  (i  to  be  large  (7)  can  be  approximated  as 

Vo  =  Vs  -  R2Iq  ~  VBE  ~  R\IC^Ri'0,Vt 


(7) 

(8) 


Eq.  (8)  is  of  the  same  type  as  that  of  the  mass  transportation 
phenomenon  described  by  (2). 

In  the  PSPICE  model  of  the  B  JT,  the  area  of  the  device  can 
be  increased  (by  300  times)  to  match  the  large  current  of  the 
fuel  cell  stack.  The  other  parameters  of  the  BJT  in  PSPICE 
library  do  not  have  much  effect  on  the  V-I  characteristic  of 
the  fuel  cell. 

In  the  simulation,  R\  and  R2  decide  the  shape  of  the  static 
characteristic  in  the  concentration  polarization  region  of  the 
fuel  cell.  In  particular,  R2  decides  the  current  I\ ,  which  corre¬ 
sponds  to  the  starting  point  on  the  concentration  polarization 


Q1  R2  lo 


Fig.  3.  Experimental  V-I  characteristic  and  simulated  characteristics  for 

different  N  for  activation  and  ohmic  regions.  Fig.  4-  Steady-state  equivalent  circuit  of  fuel  cell  module. 
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Fig.  5.  Calculation  of  R\  and  R2  from  illustrative  V-I  characteristic. 

region  of  the  V-I  curve  (Fig.  5).  The  resistance  R\  decides 
the  rate  of  change  of  voltage  in  the  concentration  polarization 
region  for  output  currents  beyond  I\ .  As  shown  in  the  illus¬ 
trative  characteristic  of  Fig.  5,  the  V-I  curve  is  approximated 
by  a  set  of  straight  lines  to  get  points  A  (A)  and  B  (A).  The 
values  of  R\  and  R2  are  calculated  as  follows. 


Fig.  6.  Comparison  of  simulated  and  experimental  characteristics. 


time  to  arrive  at  its  final  equilibrium  value  [8].  Thus,  it  is 
quite  reasonable  to  use  a  capacitor  to  model  the  capacitance 
effect  resulting  from  the  charge  double  layer.  The  capacitor 
C  is  connected  to  the  voltage  source  Vs  in  the  fuel  cell. 


3. 1.3.1.  Calculation  ofR j .  In  the  static  fuel  cell  model,  R2  is 
a  series  resistance,  which  actually  decides  the  current  limit¬ 
ing  action.  Hence,  R2  is  generally  very  small.  The  resistance 
R\  decides  the  regenerative  action  in  the  two  transistor  cir¬ 
cuit.  The  values  of  the  two  resistors  are  obtained  from  the 
experimental  data.  For  commercial  PEM  fuel  cell  modules, 
it  is  not  possible  to  reach  point  B  by  conducting  a  load  test. 
Over-current  protection  circuit  shuts  down  operation  slightly 
beyond  A.  However,  the  characteristic  can  be  extrapolated  to 
get  point  B  and  the  current  I2.  Based  on  simulation  results, 
the  following  empirical  equation  relating  I2  and  R\  is  found: 

h  =  a+~M'  <9) 

The  constants  for  the  PEM  module  tested  are  obtained  as: 
<2  =  46.8  and  £  =  74. 

3. 1.3. 2.  Calculation  of  R2.  For  a  given  value  of  R\  chosen 
from  (9),  the  empirical  relationship  between  I\  and  R2  is 
obtained  as: 

'■  =  i-  <10) 

For  the  module  tested,  the  constant  c  is  calculated  as  0.65. 

3.2.  Dynamic  model 

3.2.1.  The  charge  double  layer 

Due  to  diffusion  effects  and  the  reaction  between  the  elec¬ 
trons  (in  the  electrodes)  and  the  ions  (in  the  electrolyte),  there 
is  charge  double  layer  around  the  cathode  in  a  fuel  cell.  The 
layer  of  charge  on  or  near  the  electrode  interface  behaves 
much  like  an  electrical  capacitor.  The  result  is  that  if  the 
current  suddenly  changes,  the  operating  voltage  takes  some 


3.2.2.  Undershoot  in  stack  voltage 

The  experimental  response  of  the  commercial  fuel  cell 
module  to  a  step  change  in  load  shows  a  small  undershoot. 
After  a  load  step  to  full  power,  the  air  pump  in  the  fuel  cell 
takes  about  0.5  s  to  reach  a  higher  speed  to  provide  a  corre¬ 
spondingly  higher  airflow  rate.  In  the  mean  time,  the  output 
current  increases  very  quickly  using  up  the  air  around  the 
electrodes.  It  is  the  temporary  shortage  of  the  air  that  causes 
the  undershoot  of  the  output  voltage.  This  behaves  like  a  res¬ 
onant  circuit.  A  capacitor  C  and  an  inductor  Ls  in  series  with 
the  emitter  of  Q\  (as  shown  in  Fig.  2)  are  used  to  model  the 
undershoot  in  the  stack  output  voltage.  A  small  signal  anal¬ 
ysis  of  the  transistor  circuit  of  Fig.  2  can  be  done  to  obtain 
the  transient  performance.  It  can  be  seen  that  the  value  of 
undershoot  and  the  width  are  decided  by  the  values  of  Ls  and 
C  while  the  damping  factor  is  decided  by  the  value  of  Ls,  and 
R'2  which  is  given  by: 

*2  =  *2  Ik*  (ii) 
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Fig.  7.  Simulated  waveforms  of  fuel  cell  (C=  1  F,  Ls  =  10  mH). 
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Fig.  8.  Experimental  waveforms  of  fuel  cell. 

where  rn  is  the  base-emitter  resistance  of  each  B  JT.  The  val¬ 
ues  of  C  and  Ls  are  selected  to  match  the  transient  response 
of  the  fuel  cell.  Since  the  air  pump  delay  is  about  0.5  s,  the 
value  of  LC  used  is  large. 

4.  Results 

To  validate  the  proposed  model,  several  experiments  were 
conducted  on  the  commercial  fuel  cell  module  [  10] .  The  mod¬ 
ule  is  a  fully  integrated  system  that  produces  unregulated  dc 
power  from  a  supply  of  hydrogen  (fixed  pressure  of  72psi) 
and  air.  It  has  a  rated  power  of  1200  W.  The  output  voltage 
varies  from  43  V  dc  at  system  idle  to  about  26  V  dc  at  rated 
current  (46  A).  The  load  current  was  varied  from  a  small 
value  to  a  maximum  of  70  A  (load  resistance  in  the  range 
0.5-100  £2).  The  values  of  I\  and  h  are  62  and  77  A,  re¬ 
spectively.  The  corresponding  values  of  R\  and  R2  are  5  and 
0.01  Q.  Fig.  6  shows  the  static  characteristics  obtained  from 
simulation  and  experiment  and  the  two  agree.  For  the  dy¬ 
namic  response,  the  load  was  changed  suddenly  from  108  to 
0.55  £2.  The  simulation  and  experimental  results  are  shown 
in  Figs.  7  and  8.  It  is  seen  that  the  undershoots  in  the  two 
waveforms  agree. 

5.  Conclusions 

A  novel  circuit  model  for  the  PEM  fuel  cell  stack  is  pre¬ 
sented.  The  simple  model,  which  can  be  used  in  the  design 


of  fuel  cell  power  systems  uses  a  diode  and  a  pair  of  BJTs 
for  describing  the  static  conditions.  A  capacitor  and  an  in¬ 
ductor  are  used  to  represent  the  dynamic  conditions,  which 
occur  in  switching  power  systems.  All  the  elements  used  in 
the  model  are  from  the  PSPICE  library  and  slight  changes 
are  made  in  the  parameter  values.  The  model  is  validated  by 
comparing  the  simulation  and  experimental  results  obtained 
on  a  commercial  fuel  cell  module. 
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